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ABSTRACT 

Given a lowpass filter, finding a dual lowpass filter is an essen- 
tial step in constructing non-redundant wavelet filter banks. 
Obtaining dual lowpass filters is not an easy task. In this pa- 
per, we introduce a new method called committee algorithm 
that builds a dual filter straightforwardly from two easily- 
constructible lowpass filters. It allows to design a wide range 
of new wavelet filter banks. An example based on the family 
of Burt-Adelson's 1-D Laplacian filters is given. 

Index Terms — Laplacian pyramid, non-redundant filter 
bank, polyphase representation, wavelet filter bank 

1. INTRODUCTION 



not be biorthogonal to h. We show that given such a pair 
of lowpass filters / and g, the dual filter d that satisfies the 
above conditions (i) and (ii) simultaneously can be obtained 
straightforwardly. Since finding the pair (/, g) that satisfies 
(a) and (b) independently is much easier than finding the dual 
d directly, our method provides an easier algorithm to con- 
struct non-redundant wavelet FBs. We refer to this algorithm 
as the committee algorithm. 

In order to develop our committee algorithm we use the 
Laplacian pyramid (LP) algorithm and its polyphase repre- 
sentation. We briefly review these along with other relevant 
concepts in Section 2. Our main results are presented in Sec- 
tion 3 together with an example illustrating our findings. We 
summarize our results in Section 4. 
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A filter bank (FB) consists of the analysis bank and the syn- 
thesis bank, which are collections of, say p, filters linked by 
downsampling and upsampling operators, respectively fl] . A 
FB is typically referred to as a wavelet FB if each of its analy- 
sis and synthesis banks has exactly one lowpass filter and the 
rest of them are all highpass filters. Designing non-redundant 
wavelet FBs is an important problem since it leads to the con- 
struction of biorthogonal (or Riesz) wavelet bases under well- 
understood constraints flUH . 

Such a design problem is often reduced to the problem 
of extending a matrix with Laurent polynomial entries JfQ. In 
this approach, an essential step is to find a dual lowpass filter d 
of a given lowpass filter h so that (i) h and d are biorthogonal, 
(ii) d has positive accuracjfl (cf. see Section 2 for definition), 
and optionally (iii) d has some other properties desirable for 
the specific design. Although a filter that is only biorthogonal 
to h, or a filter that has only positive accuracy can be obtained 
readily, finding the filter that satisfies both properties (i) and 
(ii) at the same time is not easy. The filters with biorthogo- 
nality property but not with positive accuracy lead to FBs but 
not to wavelet FBs. 

In this paper we show that the dual filter d can be obtained 
by first finding two lowpass filters / and g that have a "partial 
responsibility": (a) / is biorthogonal to h but it needs not have 
positive accuracy and (b) g has positive accuracy but it needs 



'in designing non-redundant wavelet FBs, having lowpass filters with 
positive accuracy is essential as the positive accuracy guarantees that all the 
other filters are highpass (5J. 



2. LP AND ITS POLYPHASE REPRESENTATION 

We consider only the FBs with finite impulse response (FIR) 
filters that have the perfect reconstruction (PR) property. FIR 
FBs are useful as they provide fast algorithms. We use A to 
denote the n x n dilation matrix for sampling, where n is 
the spatial dimension. The PR property holds only if p > 
q := | dot A|, where p is the number of filters in each bank as 
before. If p = q, the FB is called non- redundant; otherwise, 
it is called redundant. 

The Laplacian pyramid is introduced by Burt and Adel- 
son |6). The LP algorithm has analysis and synthesis pro- 
cesses. One begins the LP analysis process with two FIR low- 
pass filters: the compression and prediction filters. Using the 
compression filter h, one obtains the coarse coefficients that 
approximate the input signal. The prediction filter g is used 
for predicting the original signal from the coarse coefficients. 
One computes the detail coefficients by subtracting the pre- 
dicted signal from the input signal, and then stores them along 
with the coarse coefficients in place of the original signal. The 
standard LP synthesis process recovers the original signal by 
reversing the above subtraction step. It is well-known that the 
LP representation is redundant; there are more coefficients 
after the LP analysis process than those in the input signal. 

The polyphase decomposition in |7| is a widely used 
method of transforming a filter (or signal) into q = | det A| 
filters (or signals) running at the sampling rate 1/q. Below 
we briefly review the polyphase representation of LP and 



refer [8][9] for details. We let T be a complete set of represen- 
tatives of the distinct cosets of the quotient group Z n /AZ™ 
containing 0, and T* be a complete set of representatives 
of the distinct cosets of 27r(((A T )~ 1 Z")/Z") containing 0. 
Here A T is used to denote the transpose of A. Then both the 
sets T and T* have q = | det A| elements. We use 

v a = 0,^i, • • • ,v q -\ 

to denote the elements of T. 

The (polyphase) LP filters are 

G(z):=[G V0 (z),G Vl (z),...,G v ,_ 1 (z)f r , 
; — [Hu (z), H Ul (z), . . . , H Uq l (z)], 

where G v and iJ„ are the z-transforms iflOl of g v and h v , 
respectively, with 

g v {m) ;= g(km+v), h v (ra) := h(Am—v), y-mGl/ 1 , 
and 



S„(z) := [ G(z) I q ] , (1) 

are the (polyphase) analysis and trivial synthesis operators of 
LP, respectively. Here and below we denote the m x m iden- 
tity matrix using I m . Since h and g are FIR filters, the entries 
of the above operators are all Laurent polynomials in z. 

In this paper we normalize the lowpass filter so that its 
coefficients sum up to Jq. The LP compression filter is called 
interpolatory if Hq(z) = 1 I Jq. We also assume that both LP 
filters have positive accuracy. We recall that for an integer 
m > 0, the filter h has accuracy m if the number of zeros of 
H(e lu) ), the Fourier transform of h, at u G r*\{0} is m (TJ. 

The pair (Alp(z), Sq(^)) of the above LP operators can be 
considered as the polyphase representation of a FB. The fact 
that LP is redundant can be seen also from these polyphase 
representations. The analysis operator Axp(z) and the trivial 
synthesis operator So (z) clearly satisfy the perfect reconstruc- 
tion property Sq(z)A.lp(z) = I q , but it is oversampled since 
h-Lp(z) and Sq(z) are not square matrices. They are of size 
(q + 1) X q and q X (q + 1), respectively. 

For a given analysis operator Alp(z), there exist infinitely 
many different synthesis operators that satisfy the perfect re- 
construction property. It is well-known ( |8]|9][IT|[T2I|) that 
the most general LP synthesis operator Slp(z) is given as the 
form 

S LP (z) := [ G(z) + V(z)B(z) I q - V(z)H(z) ] , (2) 

with B(z) := 1 - H(z)G(z), V(z) := U c (z) + U d (z)G(z), 
where U c (z) and Ud(z) are the parameter matrices of size qx 1 
and q x q, respectively, consisting of the Laurent polynomial 
entries. We note that the LP filters are biorthogonal if and 
only if B(z) = 0. 



The synthesis operator Sq(z) in (Q~|i corresponds to the 
case where V(z) = 0. The LP pair (Ai,p(z), So(z)), with 
the trivial synthesis operator, is not a wavelet FB since none 
of the filters in the synthesis bank S (z) is highpass. 

There are a few existing methods for designing wavelet 
FBs from LPs. The wavelet frame theory was used in l8l [T3l 
to reinterpret the LP process as a wavelet process, and the as- 
sociated wavelet FB was constructed along the way. However 
the resulting wavelet FB is redundant and the filters are not 
necessarily FIR unless the LP filters are biorthogonal. The 
paper [9| concerns about transforming the given LP to a criti- 
cal representation (CLP). A CLP of a given LP is said to exist 
if there exists a q x (q + 1) matrix R(z), called the reduc- 
tion operator, whose entries are the Laurent polynomials in z 
such that (i) R(z)A^p(z) is invertible, and (ii) the first row of 
R(z)Alp(z) is the same as H(z). We also recall that a CLP is 
called ECLP if the reduction operator takes a simple special 
form. The ECLP enables us to obtain a critical representa- 
tion effortlessly by eliminating a redundant portion of the LP. 
In the same article, it is shown that a special class of ECLP 
called interpolatory is available whenever the LP compression 
filter satisfies 

H vf,-ii z ) = VV?' for some /3 e (3) 

which includes the LP compression filter h being interpola- 
tory as a special case, is satisfied. It turns out that the inter- 
polatory ECLP also provides a way to design non-redundant 
FIR wavelet FBs from LPs with the condition (O. 

3. COMMITTEE ALGORITHM 
3.1. Theory 

It is easy to see that a CLP exists if and only if there exists 
a non-redundant FIR FB (Aclp(z), Sclp(z)) such that the 
first row of Aclp(z) is the same as the given LP compres- 
sion filter H(z). It is also easy to see that if such a non- 
redundant FIR FB exists, then there exists Laurent polyno- 
mials f Uo (z), f Ul (z), ■ ■ ■ , f Vq _ x (z) such that 

f V0 (z)H V0 (z) + f Vl (z)H Ul (z) + --- + f Vq _ 1 (z)H V9 _ 1 (z) = 1. 

(4) 

Since F(z) := [f VQ (z), f Vl (z), f Vg _ 1 (z)] T can be con- 
sidered as the polyphase representation of a filter, say /, the 
above condition between F(z) and H(z) is equivalent to the 
biorthgonality of the filters / and h. Furthermore it is easy to 
see that the filter / whose polyphase representation satisfies 
(HI is necessarily lowpass (see, for example, [9]). 

The next theorem shows that the above condition is also 
sufficient to obtain a wavelet FB and that there is a simple 
algorithm which produces the wavelet FB. 

Theorem 1 Suppose that the LP compression filter satisfies 
for some Laurent polynomials f„ {z), f Vl (z), ■ ■ ■ , f lJql (z). 
Then there exists a CLP with a wavelet FB. 



Proof: LetV(z) = [f Vo (z),fu 1 (z) 1 ---,f Vt _ 1 (z)] T . 
Then, from det(l, - V(z)H(z)) = 1 - H(z)V(z) = and 
^>lp{z)^lp{z) — I q , and from the form of Slp(z) given in 
(0), the rank ofI q — V(z)H(z) has tobeq — 1. Thus, there are 
exactly q — 1 linearly independent columns ofI q — V(z)H(z) 
and the other column can be written as a linear combination 
of these q — 1 columns. Suppose that the k-th column can 
be written as a linear combination of the rest of the columns. 
Then there exists a coefficient vector 



Ci,--- , Cfe_iJ 



:= [Cfc+i, 



such that (l q -V(z)K(z))[c', l,c"] T = 0. Thus the (It- 
column of 



SO) = S LP (z) 



is zero vector. Since S£p(z)Alp(z) = I g and 
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if we define Sclp{z) as the submatrix of S(z) obtained by 



deleting its (fc -+ 
of 

A(z) = 



1) -th column, and Aclp(z) as the submatrix 
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by deleting its (fc + l)-th row, we have Sclp(z)A-clp{z) = 
I q . It is easy to see that the first row of Kclp(z) is H(z). 
Hence the CLP with the FB (Pl C lp{z), Sclp(z)) exists. 

Let us now show that the FB (&clp(z),Sclp(z)) is a 
wavelet FB. Let d be the filter associated with the polyphase 
representation G(z) + V(z)B(z). Since the FB is non- 
redundant, it suffices to show that the filters h and d have 
positive accuracy (see for example [5]). Since h is assumed 
to have positive accuracy, all we need to show is that d has 
positive accuracy, which is equivalent to (cf. Result 2 in jf9|? ) 
G(l) + V(±)B(1) = ^=1 where 1 = [1, ■ • • , 1] T € W. This 
is true since, no matter what the value V(l) is, we have 



G(1) + V(1)B(1) = G(1) 

1 , w It 
:1 + V(1)(1 -V 



V(1)(1-H(1)G(1)) 
1 



where the fact that the LP filters h and g have positive accu- 
racy is used for the second equality. 

It is easy to see that if the LP compression filter satisfies 
the condition (0 for the interpolatory ECLP, then the set of 
Laurent polynomials {f v (z) : v € T} satisfying (01 always 
exists as one can choose f» ft _ 1 (z) = ^fq and f Va _ x {z) = 
for all a p. Hence our result here can be considered as a 
generalization of the interpolatory ECLP method in (9)- 



3.2. Algorithm 

Theorem 1 provides a new method to design (non-redundant 
FIR) wavelet FBs. In particular we obtain an algorithm to 
construct a dual lowpass filter from a given lowpass filter. 

Algorithm: Constructing a dual filter of a given lowpass filter 
Input: h, a lowpass filter with positive accuracy. 
Output: d, a dual lowpass filter with positive accuracy that is 
biorthogonal to h. 

Step 1: Find a lowpass filter / that is biorthogonal to h. 
Step 2: Choose a lowpass filter g with positive accuracy. 
Step 3: Compute the polyphase representation F(z), G(z) (as 
column vectors), and H(z) (as a row vector) of /, g, and h. 
Step 4: Set d as the filter whose polyphase representation is 
G(z)+F(z)(l-H(z)G(z)). 

We refer to the above as the committee algorithm as the 
filters / and g (and the input h) work together to build the dual 
of h. The filter / in Step 1 can be found, for example, by using 
the technique of Grobner bases 1 14]. There are many possible 
choices for the filter g in Step 2, and one possibility is to take 
the input h as g. Once specific / and g are chosen in Steps 1 
and 2, the dual filter d in Step 4 is uniquely determined. 

We recall that for a highpass filter I, the number of zeros 
of L(e tul ), the Fourier transform of I, at uj = is referred 
to as the number of (discrete) vanishing moments of the filter 
I (5). It is well-known (see, for example, |5 |) that for a non- 
redundant FIR wavelet FB associated with the above pair of 
lowpass filters h and d with accuracy a and j3 respectively, the 
synthesis wavelet filters have at least a vanishing moments, 
and the analysis wavelet filters have at least f3 vanishing mo- 
ments. Since d is determined by /, g and h, expressing j3 in 
terms of the properties of more attainable /, g and h is desir- 
able if possible. The next theorem shows that it can be done 
in a simple way. We omit its proof as it can be obtained by 
scrutinizing the proof of Theorem 2 in J9j- 

Theorem 2 Let Pi and P2 be positive integers. Suppose 
that the accuracy of g is f3\, and the number of zeros of 
q - H(e iM )G(e 2UJ ) at u = is p\ where H(e lu ) and G(e lu ) 
are Fourier transforms ofh and g, respectively. Then the dual 
filter d in the committee algorithm has at least min{/3i, $2} 
accuracy, hence the analysis wavelet filters in the wavelet 
FB generated by the pair (h,d) have at least min{/3i, ^2} 
vanishing moments. 

We finish this section by illustrating our findings using an 
example. 

Example: Let us consider the 1-D dyadic case with the Gaus- 
sian filters suggested by Burt and Adelson [6|. We assume 
that we are given the lowpass filter h of the form h = \/2 h 
where 

h(0) = a, h{-l) = h{\) = \, fe(-2) = fc(2) = ±-|, 



Scaling function for h with a-0.6 Scaling function for d with a=0.6 
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Fig. 1. Scaling functions associated with the filters h and d 
for a = 0.6 in the example. 

where a is a scalar parameter. For this lowpass filter, a sim- 
ple computation shows that, as long as a ^ 0.25, the Laurent 
polynomials f Q (z) = ^2 ^ and h{z) = y/2 fs=i (« + 1) 
satisfy the condition (|4j. Below we assume that a ^ 0.25. 
Although the filter / associated with the polyphase repre- 
sentation [fo(z),fi(z)] T is always lowpass, its accuracy is 
zero except the case a — 0.75, hence it cannot be used as 
the dual lowpass filter for constructing wavelet FBs directly 
unless a = 0.75. However in the committee algorithm, an- 
other lowpass filter g that is responsible for the accuracy can 
be utilized. For simplicity, we choose g = h. In this case 
mini/?!, /3 2 } in Theorem 2 is two for every a 7^ 0.25. For 
example, for the choice of a = 0.6, the dual lowpass filter d 
obtained by the committee algorithm has eleven taps. In this 
case, both h and d generate the scaling functions that are in 
L 2 (M), whose graphs are shown in Figure 1. We note that the 
choice of a = 0.5 leads to the interpolatory filter h, which is 
studied in J9). 



4. SUMMARY 

In this paper we presented a new technique called the com- 
mittee algorithm, which suggests an easy way to construct 
a dual lowpass filter of a given lowpass filter h. Finding a 
dual lowpass filter d that is biorthogonal to h and has posi- 
tive accuracy at the same time is a crucial step in constructing 
non-redundant wavelet FBs, but it is not easy. In our method, 
one needs to find two lowpass filters separately (one respon- 
sible for biorthogonality and the other for positive accuracy), 
which is much easier than finding one lowpass filter that sat- 
isfies both. Then the dual lowpass filter d can be obtained in 
a straightforward manner. 
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